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ABSTRACT 

We present an adaptation of the rotation-corrected, m-averaged spectrum technique designed to 
observe low signal-to-noise-ratio, low-frequency solar p modes. The frequency shift of each of the 
21 + 1 m spectra of a given (n, I) multiplet is chosen that maximizes the likelihood of the m-averaged 
spectrum. A high signal-to-noise ratio can result from combining individual low signal-to-noise-ratio, 
individual-m spectra, none of which would yield a strong enough peak to measure. We apply the 
technique to GONG and MDI data and show that it allows us to measure modes with lower frequencies 
than those obtained with classic peak-fitting analysis of the individual-m spectra. We measure their 
central frequencies, splittings, asymmetries, lifetimes, and amplitudes. The low-frequency, low- and 
intermediate-angular degrees rendered accessible by this new method correspond to modes that are 
sensitive to the deep solar interior down to the core {I < 3) and to the radiative interior (4 < / < 35). 
Moreover, the low-frequency modes have deeper upper turning points, and are thus less sensitive to 
the turbulence and magnetic fields of the outer layers, as well as uncertainties in the nature of the 
external boundary condition. As a result of their longer lifetimes (narrower linewidths) at the same 
signal-to-noise ratio the determination of the frequencies of lower-frequency modes is more accurate, 
and the resulting inversions should be more precise. 

Subject headings: methods: data analysis. Sun: interior. Sun: oscillations 



1. INTRODUCTION 

Our knowledge of the structure and dynamics of the 
solar interior has been considerably improved by the 
use of measurements of the properties of the normal 
modes of oscillation of the Sun. However, the Sun's 
interior is far from being fully understood, and better 
measurements of the mode parameters will also help 
to better understand the mode excitation and damp- 
ing mechanisms as well as the physical properties of 
the outer layers by better constraining the turbulence 
models. A large number of predicted acoustic oscilla- 
tion modes, defined by their radial orders (n) and their 
angular degrees (I), are not yet observed in the low- 
frequency range (i.e., approximately below 1800 /xHz), 
because the amplitude of the acoustic modes decreases 
as the mode inertia increases as the frequency decreases, 
while the solar noise from incoherent, convectivc motions 
increases: thus the signal-to-noise ratio (SNR) of those 
modes is progressively reduced. Moreover, these low- 
frequency p modes have very long lifetimes, as much as 
several years, which results in very narrow linewidths, 
hence precise frequency measurements. Thanks to 
the long-duration helioseismic observations collected by 
the spa ced-based instrument s Michelson Doppler Imager 
(MDI) (jScherrer et~al] 1X9951) and Global Oscil lations at 
Low Frequencies (GOLF) jCabriel et al.lll995D onboard 
the Solar and Heliospheric Observatory (SOHO) space- 
craft, and by the ground-based, m ulti-site Global Os- 
cillation Network Group (GONG) (|Harvev et al.lll996[ ) 
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and Birmingham S olar Oscillations Network (BiSON) 
(jChaplin et al.l[T996l ) , the frequency resolution is continu- 
ously improving and the observation of lower radial-order 
solar p modes is becoming possible. Their precise mode 
parameter determination is of great interest for improv- 
ing our resolution throughout the solar interior because 
they cover a broad range of horizontal phase velocity, and 
thus a broad range of depths of penetration. Moreover, 
these low-frequency modes have lower reflection points 
in the outer part of the Sun, which make them less sen- 
sitive to the turbulence and the magnetic fields in the 
outer layers, where the physics is poorly understood. 

The usual mode-fitting analysis consists of fitting 
the 21 -\- 1 individual-m spectra of a given multiplet 
(rijl), either individually or simultaneously. Such fit- 
ting methods fail to obtain reliable estimates of the 
mode parameters when the SNR of the individual-m 
spectra is low. Instead, various pattern-recognition 
techniques have been developed in an effort to reveal 
the presence of modes in the low-frequeiicy range (see 



e.g..lSchoulll998tlAppourchaux et al .ll2000l:ICha plin et all 



120021 : IB roomhaU et al.ll2007l and references therein). In 
the case of spatially-resolved helioseismic data (such as 
GONG and MDI observations), m-averaged spectra ap- 
peared to be a powerful tool, since for a given multiplet 
(n, I), there exist 2/ + 1 individual-m spectra, which can 
result in an average spectrum with a SNR ^ 1 once 
the individual-m spectra are corrected for the rotation- 
and structure-induced frequency shifts. The m-averaged 
spectra were employed early in the development of helio- 
seismology bv lBrownl (|1985l ). but were replaced by fitting 
the m spectra individually as the quality and the SNR of 
the data improved. However, years later, in order to fully 
take advantage of the long-duration helioseismic GONG 
and MDI instruments and reach lower frequencies i n the 
solar oscillation spectrum, iSchod (fl99aL [200l [200l and 
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TABLE 1 

Details of the long GONG and MDI analyzed time series 



Start End Fill 

Scries Date Date (%) 

3960-day GONG 1995 May 7 2006 Mar 9 84.6 
2088-day GONG 1996 May 1 2002 Jan 17 83.7 
2088-day MDI 1996 May 1 2002 Jan 17 88.9 



lAppourchaux et all ()2000f ) used the m-averaged spectra 
corrected by the modeled solar rotation to detect new low 
radial-order p modes and to set upper limits on the de- 
tcctability of the g modes. These authors demonstrated 
the potential advantage of such rotation-corrected, m- 
averaged spectra. 

We present here an adaptation of the m-averaged spec- 
trum technique in which the m-dependent shift param- 
eters are determined by maximizing the quality of the 
resulting average spectrum. The analysis is performed 
on long-duration time series of the spatially-resolved he- 
lioseismic GONG and MDI observations of the low- and 
medium-angular degrees (1 < / < 35). This range of os- 
cillation multiplets samples the radiative interior down 
to the solar core. In Sec. [21 we introduce the different 
datasets used in this analysis. In Sec.[3l we describe this 
new technique in order to observe low signal-to-noise- 
ratio, low-frequency p modes, explaining the different 
steps of the analysis from the mode detection to peak- 
fitting. In Sec. m we demonstrate that this method al- 
lows us to successfully measure lower-frequency modes 
than those obtained from classic peak-fitting analysis of 
the individual- m spectra by comparing with other mea- 
surements obtained from coeval datasets. In Sec. [HI we 
present the mode parameters of these long-lived, low- 
frequency acoustic modes down to w 850 fiHz extracted 
from the analysis of 3960 days of GONG observations us- 
ing the m-averaged spectrum technique. Finally, Sec. [6] 
summarizes our conclusions. 

2. OBSERVATIONS 

Details of the spatially-resolved helioseismic observa- 
tions collected by both GONG and MDI used for this 
work (the starting and ending dates, and their corre- 
sponding duty cycles) are given in Table [TJ Coeval 2088- 
day observations of GONG and MDI were analyzed for 
oscillation multiplets with angular degrees from / — 1 
to I = 35, and a re then directly compared to those of 
iKorzenniM (|2005f l for I < 25 measurements of the same 
datasets. We also applied the analysis to 3960 days of 
GONG data (1 < / < 35), which constitutes so far the 
longest time series (w 11 years, spanning most of solar 
cycle 23) of spatially-resolved observations analyzed. 

3. METHOD 

An m-averaged spectrum corresponds to the average of 
the 21 + 1 individual-m components of an oscillation mul- 
tiplet {n,l), thus reducing the non-coherent noise. Be- 
fore averaging, each m spectrum of a given mode (n, I) 
is shifted by a frequency that compensates for the effect 
of differential rotation and structural effects on the fre- 
quencies. The m-averaged spectrum concentrates, for a 
given multiplct (n, I), all of the 2/-I-1 m components, as it 



would be if the Sun were a purely-spherical, non-rotating 
object. Thus, the average of the 2/ + 1 individual-m 
spectra considerably improves the SNR of the resulting 
m-averaged spectrum. 

3.1. Determination of the shifts 

The m-averaged spectrum is obtained by finding the 
estimates of the splitting coefficients, commonly called 
a-coefficients, which maximizes the likelihood of the m- 
averaged spectrum. The a-coefficients are individually 
estimated through an iterative process, with the initial 
values taken from a model. Thus, for a given mode (n, /), 
the frequency shift Svnim is parameterized by a set of 
coefficients, as: 

i=l 

where ai{n, I) are the splitting coefficients, and P/„j cor- 
responds t o the Clebsch-Gordan polyn omial expansion as 
defined bv lRitzwoUer fc Lavelvl (|l991f ). In this definition, 
the odd orders of the a-coefRcients describe the effects of 
solar rotation, while the even orders correspond to depar- 
tures from spherical symmetry in the solar structure as 
well as to quadratic effects of rotation. Each is chosen 
to maximize the likelihood of the m-averaged spectrum. 
This is performed through an iterative procedure. For a 
particular order i of the coefficients a,;, a range of values 
is scanned around its initial value, while the other a^/^iS 
are kept fixed to their previously estimated values. 

For each scanned value of , the individual-m spectra 
are shifted by the corresponding Clebsch-Gordan poly- 
nomials, and the mean of these 21 + 1 shifted spectra is 
taken. The mean power spectrum is then fitted using a 
Maximum-Likelihood Estimator (MLE) minimization as 
described in Sec 13.21 and its likelihood determined. 

For a Monte-Carlo simulation, the left panel of Fig. [T] 
shows the variation of the likelihood from the MLE 
minimization as a function of the first splitting coeffi- 
cient fli, showing a well defined minimum which repre- 
sents the best value of ai. The artificial power spectra 
were simulated following t he methodology described in 
iFierrv Fraillon et al.l ()1998l ). We also examined the sen- 
sitivity of the mode linewidth and the entropy as criteria 
for deterinining the best shifts. In our case, the entropy 
(|Shannonl Il948f ) can be seen as a measure of random- 
ness in the m-averaged spectrum, S, and is defined as 
-El^xln^]. 

Both linewidth and entropy show well defined minima 
around the input value of ai (middle and right panels 
of Fig. [1] respectively). Indeed, the m-averaged spec- 
trum gets narrower as ai converges to its input values 
of 0.4 ^Hz and ai = 400 nHz. Similar variations are ob- 
tained for all the a^s. As detailed in Appendix 1X1 the use 
of these different criteria to determine the best estimates 
of the a-coeficients returns consistent results. 

The iteration is performed until the difference between 
two iterations in each of the computed ai coefficients falls 
below a given threshold (such as 0.25(t in the case of ai). 
Also, in order to remove any outliers, some quality checks 
are performed after each measure of an which needs to 
fall within a constrained range of values. For example, a 
±15% window around its theoretical expectation is used 
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Fig. 1. — — log(likelihood) (left panel), natural logarithm of the mode liiiewidth {middle panel), and entropy {right panel) of the m- 
averaged spectrum as a function of the shift coefficient ai obtained from ten simulated / = 15 spectra. The vertical dashed lines represented 
the introduced value of ai (400 nHz) in the simulations. The horizontal dotted-dashcd line in the center panel represents the mode lincwidth 
of the simulated spectra, 0.4 /^Hz. 



for fli . Here, we fitted only the six first Ui in the Clebsch- 
Gordan expansion, even though the quality of the data 
supports the determination of higher-order coefficients. 

Finally, low SNR peaks in the m-averaged spectrum 
(after adjustment) are tested against the HO hypothesis. 
In the framework of that hypothesis, the resulting spec- 
tra are tested against a statistics pertaining to pure noise 
(x^ with 2(2Z-|-1) d.o.f). This test has been widely ap- 
plied to helioseismic observations in the search for long- 
lived, low radial-order p modes and g modes (see e.g., 
I Appourchaux et al] 120001 ) . In the present analysis, we 
rejected peaks that have a greater than 10% chance of 
being due to noise in the 238 analyzed windows, each 
containing 288 frequency bins. Here the fixed number of 
bins was chosen because we know that the range of theo- 
retical frequency lie within 1.5 /iHz or so. Figures [2] and 
[3] illustrate the advantage of using the m-averagcd spec- 
trum technique in the case of two oscillation multiplets 
for 2088 days of GONG data, where the m-averaged spec- 
tra before and after the correction for the splitting coeffi- 
cients are shown. These examples show the m-averaged 
spectra of the modes Z = 3, n = 5 at w 1015.0 /iHz 
(Fig. El, and Z = 16, 71 = 4 at « 1293.8 /iHz (Fig. El), 
as well as the corresponding m — v diagrams. These two 
examples were chosen to demonstrate the performance 
at different SNR levels. The corresponding 10% proba- 
bility levels are given. The m — v diagrams in the case 
of the mode Z = 3, n = 5 (right panels on Fig. ^ do 
not show any high SNR structure before or after correc- 
tion. However, the m-avcraged spectrum after correc- 
tion clearly shows the target mode (lower left-panel on 
Fig- El) I with an unambiguous detection level. The mode 
/ = 16, n = 4 presents a higher SNR (Fig. |3l) and its 
m — V diagram shows that the individual-m spectra line 
up after correction (lower-right panel on Fig.EJ. The es- 
timated splitting coefficients of the low-frequency modes 
with 1 < Z < 35 measured in the 3960-day GONG dataset 
are shown in Fig. [4] as a function of frequency and v/L 
(with L = \Jl{l + 1)), which is approximately propor- 



tional to the sound speed at the mode's inner turning 
point. Modes with selected ranges of radial orders are 
represented with different colors and symbols. 

3.2. Extraction of the mode parameters 

For a given mode {n,l), the best estimates of the 
splitting coefficients determined as discussed in Sec 13.11 
are used to calculate its m-avcraged spectrum. When 
N independent power spectra are averaged together, 
the statistics of the mean power spectrum corresponds 
to a with 2 x i V deg ree-of-freedom (d.o.f.) statis- 
tics. lAppourchauxl (|2003[ ) demonstrated that the mean 

01 21 + \ independent power densities, which has a 
with more than 2 d.o.f. statistics, can be correctly fit- 
ted with a Maximum-Likelihood Estimator (MLE) min- 
imization code developed for spectra following a with 

2 d.o.f. statistics. The asyn imetric Lorentzian model 
of iNigam fc KosovichevI (|1998[ ) was used to describe the 
m-averaged spectrum, as: 



(1 + an,lXn,lY + < 



where 



n,l' 



^jv ~ Vn,l) 



(2) 



(3) 



Then, for a given mode (n, I), the central frequency, the 
FuU-Width-at-Half-Maximum (fwhm), and the power 
height of the spectral density are respectively Vn,u ^n,u 
and Hn^i- The peak asymmetry is described by the pa- 
rameter a„.i, while represents an additive, constant 
background level in the fitted window. The first spatial 
leaks {51 = 0, 5m = ±2), commonly called m-leaks, are 
also included in the fitting model and added to Eq. El 
The frequencies of the m-leaks are set from the central 
frequency of the target mode using the previously mea- 
sured splitting coefficients (Sec. 13. ill . Their peak asym- 
metries are assumed to be the same as that of the target 
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Fig. 2. — Examples of m-avcragcd power spectra (left panels) before (top) and after (bottom) correcting for the shift coefficients in the 
case of the oscillation multiplet Z = 3, n = 5 observed in 2088 days of GONG data. The corresponding m — u diagrams are also shown 
{right pa nels). The crosses indicate the posit ion of the corresponding theoretical central frequency calculated from Christensen-Dalsgaard's 
model S I IChristensen-Dalsgaard et al.lll996l1 . The dot-dashed lines on the left hand-side panels give the 10% probability limit that a peak 
is due to noise in the 238 windows, 1.5 /xHz wide. The illustrated spectral window contains the 2/ -f- 1 components of the represented 
multiplet. 
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Fig. 3. — Examples of m-averagcd power spectra (left panels) before (top) and after (bottom) correcting for the shift coefficients in the 
case of the oscillation multiplet / = 16, n = 4 observed in 2088 days of GONG data. The corresponding m — v diagrams are also shown 
(right pa nels). The crosses indicate the posit ion of the corresponding theoretical central frequency calculated from Christensen-Dalsgaard's 
model S l IChristensen-Dalsgaard et al.lll996IV The dot-dashed lines on the left hand-side panels give the 10% probability limit that a peak 
is due to noise in the 238 windows, 1.5 /iHz wide. The illustrated spectral window contains the 2/ + 1 components of the represented 
multiplet. 



mode, while their FWHMs are a free parameter of the fit 
and different from the target mode. The amplitude of 
the m-leaks is specified to be a fixed fraction of the cen- 
tral peak, which was estimated from the leakage matrix 
developed especially for the GONG (jHiU fc Howdfl99l 
and MDI (J. Schou, private communication) data. The 
first spatial leaks in the m-averaged spectrum were deter- 
mined by averaging for a given multiplet (n, I) the (5m ±2 
leaks over the entire 2Z + 1 spectra. 



The size of the fitting window, fi^, is proportional to 
the first estimates of the mode width, r„_;, and centered 
around the frequency of the target mode. It is defined 
as: 

r!, = 20y^r2 ;-f Ai/2-^A5™, (4) 

where l^Vj. is the frequency resolution of the power spec- 
trum. The first spatial leaks are always included in 
the fitting range by adding the offset A^m = 800 nHz. 
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Fig. 4. — First six splitting coefficients of the low- frequency p modes with 1 < I < 35 measured in 3960 days of GONG data. The 
a-coefHcients are represented as a function of the mode frequency u {left column), and u/L, with L = \/l{l + 1) (right column). The 
different colors and symbols correspond to selected ranges of radial orders n: green triangles, modes with n = 1,2; purple stars, n = 3,4; 
blue upside-down triangles, n = 5,6; orange squares, n = 7, 8; and black dots, n > 9. 



The multiplicative factor 20 ensures a good sampling of 
the mode profile in the low-frequency range. A com- 
parable de f inition of the fitting window was adopted by 
iKorzennikI (|2005D . Bad fits were removed based on a set 
of quality criteria based on the fitted mode parameters 
and associated uncertainties, such as, (1) the error of the 
mode frequency must be less than its mode width; (2) 



the SNR must be larger than 1; and (3) the mode width 
must be larger than the frequency resolution. A dis- 
cussion on the impact of the fitting model (asymmetry, 
spatial leaks) on the extracted mode parameters used to 
describe the m-averaged spectrum can also be found in 
the Appendix IB] 
Figurc[5]shows examples of the m-avcraged power spec- 
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tra for four different radial orders n of the multiplct 
I = 17, and the corresponding best MLE fits, which in- 
cluded the mode asymmetry and the Sm±2 spatial leaks. 
The blending of the first m leaks is particularly clear as 
the linewidths increase with increasing frequency. 

3.3. Mode parameter and a-coefficient uncertainties 

The mode parameter uncertainties are established in 
the usual manner by the inverse of the covariance ma- 
trix. However, because the m-averaged spectrum is fit- 
ted using a MLE m inimization and, as explained in 
lApDOurchauxl (j2003[ ). the formal uncertainties must be 
normalized by the square root of the number of aver- 
aged spectra, i.e., in our case, by \/2l + 1. But this a- 
posteriori error normalization is correct only if the 2/4-1 
spectra of a given (n, I) mode have the same variance (or 
SNR). Since the condition of equal SNR among the m 
spectra within a multiplct is not satisfied in our case, the 
uncertainties of the mode parameters have to be taken as 
a first approximation only. However, Monte-Carlo simu- 
lations show that this error normalization holds even in 
the case of TO-dependent SNR (see Sec. 13.3.2)1 . 

It can also be derived that the errors on the a- 
coefhcients can be estimated as follows: 

'Ta- = ^(EK™(WO]')^.-f, (5) 

m 

where i is the a-coefficient order and P;\„ the associated 
Clcbsch-Gordan polynomials. The derivation of Eq. [5] is 
detailed in Appendix ICl 

3.3.1. m/l dependence of the signal-to-noise ratio 

Figure [6] shows the dependence in m/l of the SNR in 
the GONG data. This was obtained with modes ob- 
served in the 3960-day GONG dataset below 2000 ^Hz 
and of angular degree up to / = 35. Both mode amplitude 
and background noise depends on the azimuthal order m 
and can be described with polynomials with only even 
terms, the polynomials being different for both param- 
eters. Note that any frequency dependence of the m/l 
dependence is averaged out in Fig. [HI 

The m/l dependence of the SNR implies that the a- 
coefhcients are not exactly orthogonal and that their er- 
rors are correlated (see Appendix [C|) . However, as a first 
order approximation, the errors on the a-coefhcients can 
be estimated by using Eq. [5] (see Sec. I3.3.2p . 

3.3.2. Validation of the error estimates: Monte-Carlo 
simulations 

The formal uncertainties of the mode parame- 
ters and of the a-coefficients were verified through 
Monte-Carlo simulations. The artificial power spectra 
were simulated following t he methodology described in 
iFierrv Fraillon et al.1 () 19981 ). In a first scries of simulated 
spectra, the m dependence in amplitude within a given 
multiplct (n, /) was introduced, the SNR being symmet- 
ric in |m| around the m = spectrum. In a second series, 
no m dependence was introduced, i.e., a constant SNR 
over m. The mean values of the formal errors returned 
by the MLE minimization were compared to the RMS 
value of the corresponding fitted parameter. The Monte- 
Carlo simulations showed that in both cases the formal 



uncertainties of the m-averaged spectra determined as 
in Sec 13.31 using a MLE minimization are a very good 
approximation of the errors. 

4. COMPARISON WITH OTHER MEASUREMENTS 

4.1. Comparison with spatially-resolved 
observations (I < 25) 

GONG and MDI use two independent peak-finding ap- 
proaches to extract the mode parameters. Developed in 
the early 1990s, and mostly unchanged since, they pro- 
vide mode parameters on a routine basi s. Time series of 
108 d ays are used by the GONG project () Anderson et al.l 
|1990D. while the MDI pr oject uses 72-day time series 
(jSchoul |1999[ ). Recently, iKorzennikI (j2005[) developed 
a new and independent peak-finding method of the 
individual-m spectra, optimized to take advantage of the 
long, spatially-resolved, helioseismic time series available 
today from both projects. 

4.1.1. Mode detection: I — v diagram 

lKorzenn"i3 ()2005f ) applied his peak-fitting to extract the 
low- and medium-degree (/ < 25) mode parameters from 
both GONG and MDI observations using one 2088-day 
long time series, as well as using five overlapping seg- 
ments of 728 days. In order to compare our results ob- 
tained with the m-averaged spectrum technique, we ap- 
plied the procedure described in Sec. [3] to the same 2088 
days of GONG and MDI observations (Table [J). Fig- 
ure[7]shows the l—v diagrams of the low-frequency modes 
measured with the two different analyses in the case of 
the 2088-day GONG {left panel) and MDI (right p anel) 
datasets. The modes measured bv IKorzennikI (j2005l) with 
a classic peak-fitting method of the individual-m spectra 
are represented by the o pen circles. We co nsidered that 
a given mode (n, /) from IKorzennikI ()2005t ) was detected 
when at least two of the 21 + 1 m spectra were success- 
fully fitted, which is enough to obtain estimates of the 
corresponding central frequency and first splitting coef- 
ficient oi. The red dots represent modes measured with 
the m-av eraged spec t rum t echnique which were not ob- 
served bv IKorzennikI ()2005f ). A significantly larger num- 
ber of low-frequency modes in the 2088-day GONG and 
MDI datasets (respectively 45 and 14 new modes) down 
to « 900 /iHz can be measured using the m-averaged 
spectrum technique. 

4.1.2. Mode parameter and uncertainty comparisons 

In order to check the accuracy of the technique and 
to identify any potential bias in our analysis, we com- 
pare the central frequencies and splitting coefficients ob- 
ta ined by the t wo rn ethods. The individual-m frequencies 
of lKorzennid (pOOSi ) w ere fitted using a Clebsch-G ordan 
polynomial expansion (jRitzwoUer fc Lavelvlll99il ) in or- 
der to estimate the corresponding central frequencies and 
a-coefficicnts of each (n, /) multiplct. The formal un- 
certainties of the individual-7Ti frequencies were used as 
fitting weights. The left panel on Fig. |8| shows the dis- 
tribution of the differences in central frequencies below 
K, 1800 //Hz of the common modes between the 2088- 
day GONG estimates measu red using the m- averagcd 
spectrum technique and from IKorzennikI (|2005( l (as rep- 
resented on Fig. [7]), demonstrating that there is no 
frequency dependence over the analyzed low-frequency 
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Fig. 6. — Signal-to-noisc ratio for modes up to i = 35 as 
function of m/l, obtained with the 3960-day GONG datasct. 



range. The distribution was fitted by a Gaussian func- 
tion, and its associated parameters (mean, standard de- 
viation) are indicated on Fig. [51 While, on average, 
the GONG central frequencies obtained using the m- 
aver aged spectrum tec hnique are less than 1 nHz smaller 
than iKorzennikI ()2005f )'s estimates, this offset is not sig- 
nificant — the corresponding standard deviation being 
about 5 times larger. The MDI frequencies estimated 
with the m-avcraged spectrum technique give compara- 
ble, insignificant mean differences with Korzemiik ( 200^. 
Similar results are obtained with the splitting coeffi- 
cients. 

We also compared the low frequencies {v < ISOO^Hz, 
see Fig. [T]) estimated in both the 2088-day GONG and 
MDI datascts using the m-averagcd spectrum technique. 
The right panel of Fig. [8] represents the distribution of 
the frequency differences of the common modes, in the 
sense GONG minus MDI. The mean difference value is 
of -0.17 ± 1.99 nHz, i.e. the GONG and MDI low- 
frequency m odes are essentia lly the same. The mean 
difference in lKorzennii3 (|2005l ) 's central frequencies be- 
tween the 2088-day GONG and MDI datasets for modes 
below 1800 fiUz is of 0.35 ± 5.40 nHz. The splitting 



coefficient estimates are also consistent between the two 
datasets with in the case of the ai coefficient a mean 
difference of -0.04 ± 0.31 nHz. 

The left panel of Fig. [5] shows the formal la uncertain- 
ties of the central frequencies of the measured common 
modes between the 2088-day GONG dataset analyzed 
in the present analysis a nd th e coeval 2088-day GONG 
dataset from ]Korzcnnik^ (|2005[ ) up to « 1800 /iHz. Our 
estimates of the freq uency uncertaintie s are much smaller 
th an those quoted b v lKorzenni ^ (|2005f ). However, Fig. 10 
in IKorzennikI ()2005[ ) suggests that the errors are overesti- 
mated, and that his results "might be too conservative" . 
Korzemiik (2008,) reported that in the case of a 2088-day 
long time series, as a first estimate, a multiplicative factor 
of 0.75 need s to be apphed to t he frequency uncertainties 
reported in iKorzenniM (|2005f ). Ho wever, despite the se 
uncertainty scaling issues, while the iKorzeniiig ()2005l )'s 
uncertainties show an increase with decreasing frequency 
from sa 1500 /iHz, the uncertainties returned from the m- 
avcragcd spectrum technique do not show this increase, 
thanks to the higher SNR of the m-averaged spectrum 
than for the individual-m spectra. 

The uncertainties on the a-coefficients returned by the 
m-averaged spectrum technique are also smaller than 
the ones obtained by fitting the individual-m spectra, 
as shown on the right panel of Fig. [9] in the case of the 
ai coefficients. As for the frequencies, the a-coefficients 
of the modes below w 1500 /iHz are better constrained 
using the m-averaged spectrum technique. 

4.2. Comparison with Sun-as-a-star 
observations (I < 3) 

The spatially-resolved GONG and MDI instruments 
are not optimized to observe low-degree solar p modes 
below / < 3, unlike the Sun-as-a-star, integrated-light 
instruments such as the space-based instrument GOLF 
onboard SOHO and the ground-based, multi-site BiSON 
network. The low-degree modes are of particular interest 
as they reach the very deep interior of the Sun. How- 
ever, the spatially-resolved observations are still able to 
observe such low-degree oscillations. 

Low-degree (Z < 3) modes down to sa 1000 /iHz are 
observed in both GONG and MDI data with the m- 
averagcd spectrum technique as illustrated in Figs.[7]and 
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Fig. 8. — Left panel - Histogram of the differences (in nHz) in the estim ated centra l frequ encies between the 2088-day GONG dataset 
using the m-averaged spectrum technique and the coveal 2088-day GONG iKorzenmkj 1 120051 ) 's estimates. Right panel - Histogram of the 
differences (in nHz) in the estimated frequencies between the 2088-day GONG and MDI datasets using the m-averaged spectrum technique 
in the sense GONG minus MDI. The corresponding Gaussian function fits and their associated mean values and standard deviations are 
also indicated. 



[T2l In order to test the capability and the precision 
of the m-averaged spectrum technique to observe low- 
degree, low-frequency modes in spatially-resolved data, 
measurements obtained for w 11 years of the Sun-as- a- 
star GOLF and BiSON instruments were compared with 
the 3960-day GONG dataset and the 2088-day GONG 
and MDI datasets. The GOLF data were indepen- 
dently analyzed by two mode-fitting algorithms (R. A. 
Garcia, private communication"*, and P. Boumicr, pri- 
vate communication^). The BiSON observations come 
from a combination of frequencies obtained with dif- 
ferent long time series analyzed in order to measure 
low-frequency moc les (W. J. Chaplin, private commu- 
nicat ion. See also lChaplin et a"ni2002t iBroomhall et al.l 
|2007|) . However, while imaging instruments give us 
access to all of the 21 ~\- \ individual-m components, 
only / + 1 components can be clearly observed with 
integratcd-sunlight technique. These various compo- 
nents have different spatial structure over the solar sur- 

4 From 1996 April 11 to 2006 April 18. 

5 From 1996 April 11 to 2006 May 23. 



face, which can lead to differences in the extract e d cen - 
tral frequencies of the multiplet. [Chaplin et al.l (j2004l ) 
and lAppourchaux fc Chaplini (j2007t ) derived expressions 
to predict the differences between the low-degree fre- 
quencies extracted from spatially-resolved (as MDI and 
GONG) and Sun-as-a-star (as GOLF and BiSON) obser- 
vations. However, in the following, as a first approxima- 
tion, we compared directly the extracted mode parame- 
ters. 

The comparisons of the estimated mode frequencies 
and oi rotational splittings between the common low- 
degree (1 < ^ < 3), low-frequency modes in the two types 
of observation are shown on the left and right panels re- 
spectively of Fig. [ini The three different datasets and 
analysis methods give consistent results, for both the fre- 
quency and the splitting coefhcient oi. Of course, this is 
only assuming that the different subsets of observed mul- 
tiplets from both types of observational technique "see" 
the same central frequencies. 

Thanks to decade- long available datasets, the low- 
degree, low-frequency modes arc today measured lower 
than 1200 /-tHz with high precision, demonstrated by 
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frequency differences were normalized by the combined errors. Right panel - Same as on the left panel, but for the splitting coefficient ai . 



the consistency in the extracted parameters from differ- 
ent instruments using distinct and independent analy- 
sis. Figure [10] also demonstrates that spatially-resolved 
observations can provide as accurate measurements of 
the low-degree modes as the Sun-as-a-star instruments 
do. Moreover, the m-averaged spectrum technique allows 
the observation of lower radial-order I = 3 modes than 
the integrated-light GOLF and BiSON observations, for 
commensurate observation lengths, thanks to the obser- 
vations of the 21 + 1 components (Fig. [TT|l . 

5. MODE PARAMETERS OF THE LOW-FREQUENCY 
OSCILLATIONS 

The m-averaged spectrum technique has been applied 
to 3960 days of GONG observations (see Sec. span- 
ning most of the 11 years of solar cycle 23. The anal- 
ysis covered low-frequency modes with angular degrees 
from Z = 1 to / = 35. Oscillation multiplets well below 
1000 /iHz were detected with good precision, such as the 
modes / = 4, n = 4 at w 913.5 fiHz; I = 9, n = 3 at 
« 930.5 /iHz; / = 16, n = 2 at « 912.1 /iHz; or / = 31, 
n = 1 at ~ 907.5 /iHz. Some examples are illustrated in 
Figs. [2] and [3l These low horizontal-phase- velocity modes 



do not penetrate deeply into the Sun, but their very high 
inertias afford higher precision frequencies for the inver- 
sions. It is clear from Sec.|4]that this method allows us to 
observe modes that are otherwise lost in the background 
of each individual-m spectrum of a given multiplet (n, /), 
and thus unobservable with a classic peak-fitting anal- 
ysis. The I — V diagram of the observed low-frequency 
modes (1 < ' < 35) down to n = 1 and w 850 /iHz in 
the 3960-day GONG dataset and 2088-day GONG and 
MDI datasets with the m-averaged spectrum technique 
is shown on Fig. 

5.1. Mode linewidths, heights, and background levels 

Figure [T51 shows the fitted mode fwhms r„.; {upper- 
left panel) and mode heights Hn^i {upper-right panel) of 
the measured low-frequency oscillations. The fitted back- 
ground level is also represented on the right panel. The 
FWHMS and heights are extremely valuable tests of mod- 
els of the physical processes responsible for the mode 
damping and excitation by the turbulent convective mo- 
tions in the outer layers of the Sun: the mode damping is 
inversely related to the fwhm of the mode, and the mode 
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excitation is proportional to the mode height times the 
mode FWHM squared (for a detai led description, see, e.g., 
iSalabert fc Jimenez-Revesll2006[ ). The leveling off of the 
mode widths observed below k, 1100 /xHz, despite their 
dispersion becoming larger, could be a resolution effect, 
the peaks being then so narrow that the limiting resolu- 
tion o f the spectrum becomes an issue. Moreover. ISchod 
( |2004f ) did not observe such behaviour at low frequency 
in MDI data with a 2952-day time series. 

As indicated by different colors and symbols on Fig.fTSl 
the fitted mode widths follow ridges for equal radial or- 
ders n. This dependence on angular degree [V) is directly 
related to the mode inertia (/) in terms of a power law, 
as illustrated on the lower-left panel of Fig. [131 The I 
dependence in the mode fwhms is removed when repre- 
sented as a function of the mode inertia /. 

5.2. Mode asymmetry 

The mode parameters extracted through the routine 
GONG and MDI peak-fitting pipelines are obtained by 
use o f sym metric Lorenz tian profiles () Anderson et al.l 
(|1990( ) and ISchoul (|1992l ) respectively). However, it 
was demonstrated that ignoring the peak asymmetry in 
the description of the acoustic modes leads to biais in 
the estimated mo de parameters (see Appendix [B] and 
iThierv et al.ll200(3l ). Today, most of the estimates of the 
mode asymmetries have been restricted to low degrees 
{I < 3) only, from S un-as-a-star , integ rated-sunlight ob- 
servations. However, iKorzenniS (j2005[ ) used asymmetric 
profiles and presented estimates of the peak asymme- 
try for modes with angular degrees 1 < Z < 25, ob- 
tained with GONG a nd MDI observations. Recently, 
iLarson fc Schoul (|2008f ) are planning to reprocess all the 
MDI medium-^ data including a set of corrections and im- 
provements (such as the mode asymmetry) in the MDI 
pipeline algorithm itself. 

The asymmetry parameter {an,i) in the low-frequency 
range, obtained by fitting the 3960-day GONG m- 
averaged spectrum (1 < / < 35), is shown on the lower- 
right panel of Fig. [T21 The extracted peak asymme- 
try is well constrained down to w 1400 /iHz, with a 



mean value of about —0.044 ± 0.002, and no disccrn- 
able / dependence. The average asymmetry observed 
in the m-averaged spectrum is consistent with other 
me asurements . For instance, the mean value observed 
bv iKorzennig (j2005h was about —0.04 for modes below 
2000 /xHz and / < 25, once his estin i ates a re transformed 
back into the lNigam fc KosovichevI ()1998f )'s definition of 
the peak asymmetry. A comparable mean value is also 
observed at the lowest frequencies for which asymme- 
tries were reporte d in Sun-as-a-star, integrated-sunlight 
observations (e.g.. iThierv et al.]|200Cl[ ). 

5.3. Mode frequencies 

Figure [14] shows the frequency differences (in iiHz) be- 
tween the fitted low-frequency modes observed in the 
3960-day GONG dataset using the m-averaged spec- 
trum technique and the corresponding theoretical val- 
ues calculated from Christensen-D alsgaard's model S 
(jChristensen-Dalsgaard et al]|1996[ ). The corresponding 
frequency uncertainties were multiplied by 20 to render 
them visible. These comparisons are represented as a 
function of the angular degree {left panel), of the fre- 
quency (middle panel), and of the inner turning point 
(right panel). Modes of equal radial orders are connected. 
As these differences between observed and theoretical fre- 
quencies show, there is still room to improve the model 
of solar internal structure. Note that the right panel on 
Fig. [Til illustrates also the wide range of depths of pene- 
tration that these low-frequency modes cover. 

6. CONCLUSION AND DISCUSSION 

We presented here an adaptation of the rotation- 
corrected, m-averaged spectrum technique to observe 
low signal-to-noise-ratio, low-frequency solar p modes in 
spatially-resolved helioseismic data. For a given multi- 
plet {n,l), the shift coefficients describing the differen- 
tial rotation- and structural-induced effects are chosen 
to maximize the likelihood of the m-averaged spectra. 
The average of the 2Z -I- 1 individual-m spectra can re- 
sult in a high signal-to-noise ratio when the individual-m 
spectra have a too low signal-to-noise ratio to be suc- 
cessfully fitted. This technique was applied to long time 
series of the spatially-revolved GONG and MDI observa- 
tions for low- frequency modes (i.e., approximately below 
1800 /iHz) with low- and intermediate-angular degrees 
(1 < ^ < 35). We demonstrated that it allows us to mea- 
sure lower frequency modes than with classic peak-fitting 
analysis of the individual-m spectra. Figure [TSl shows the 
new low-frequency solar p modes observed in spatially- 
resolved data using the m-averaged spectrum technique 
in long time series of both GONG and MDI observations. 
Their central frequencies and splitting ai coefficients, as 
well as their associated uncertainties are indicated in Ta- 
ble [21 These normal modes of oscillation were predicted 
but were not measured previously. The potential of the 
m-averaged spectrum technique for increasing our knowl- 
edge of the solar interior is clearly illustrated on Fig. \TE\ 
p modes well below 1000 /xHz being measured with a 
high accuracy thanks to their longer lifetimes. We also 
demonstrated that the m-avcragcd spectrum technique 
returns unbiased results with no systematic differences 
with other long-duration measurements, which also in- 
clude the asymmetry in the mode profile description. 
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Fig. 12. — I — u diagram of the low-frequency solar p modes from i = 1 up to / = 35 observed with the m-averaged spectrum technique 
in 3960 days of GONG observations (large black circles), and 2088 coeval days of GONG (medium green circles) and MDI (red dots) 
observations. The ridges of same radial order are also indicated from n = 1 to n = 12. 
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Fig. 13. — Mode fwhm F„ j (upper-left panel) and height H„ ; (upper-right panel) as a function of frequency for the low-frequency 
p modes I < 35 extracted from 3960 days of GONG data with the m-averaged spectrum technique. The different colors and symbols 
correspond to selected ranges of radial orders n: green triangles, modes with n = 1,2; purple stars, n = 3,4; blue upside-down triangles, 
n = 5, 6; orange squares, n = 7, 8; and black dots, n > 9. The frequency resolution of the analyzed dataset is indicated on the upper-left 
hand-side plot, and the fitted background level (i?n,i) is also represented on the upper-right hand-side plot (small black dots). Lower-left 
panel - Same as above, but for the mode FWHM F„,j as a function of the mode inertia (/), calculated from Christensen-Daslgaard model S 
ijChristensen-Dalsgaard et al.|[T996I V The frequency resolution is also indicated. Lower-right panel - Same as above, but for the mode 
asymmetry (q:„^;). 
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14. — Differences between the low-frequency central frequencies observed in 3960 days of GONG data (1 < ' < 

the m-averaged spectru m technique and the associated theoretical values, taken from Christcnsen-Dalsgaard's model S 

IIChristensen-Dalsgaard et al.M1996l) . as a function of the angular degree {left panel), of the frequency {middle panel), and of the inner 
turning point {right panel). The uncertainties multiplied by a factor 20 are also presented. Ridges of same radial orders are connected, and 
for guidance some of the corresponding radial orders are indicated in the middle panel. 
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TABLE 2 

Set of the new low-frequency solar p modes observed in 
the gong and mdi datasets with the m-averaged 
spectrum technique in the range 1 < z < 35. 



I n 


Frequency 


ai-cocfficicnt 




(mHz) 


(nHz) 



1 


7 


1185.599 


± 


0.005 


431.491 




6.161 


2 


7 


1250.555 


± 


0.003 


428.263 




2.286 


3 


5 


1015.046 


± 


0.005 


430.154 


± 


2.471 


4 


4 


913.477 


± 


0.004 


420.055 




1.594 


4 


5 


1062.140 


± 


0.002 


429.275 




0.614 


5 


4 


954.560 


± 


0.002 


430.712 




0.596 


6 


4 


992.412 


± 


0.002 


431.906 




0.502 


6 


5 


1145.074 


± 


0.002 


432.019 


± 


0.464 


7 


4 


1028.156 


± 


0.003 


430.292 


± 


0.740 


8 


4 


1062.338 


± 


0.002 


434.087 




0.459 


9 


3 


930.540 


± 


0.002 


430.045 


± 


0.363 


11 


3 


987.206 


± 


0.002 


436.639 


± 


0.344 


12 


3 


1013.572 


± 


0.001 


435.158 


± 


0.172 


13 


3 


1038.795 


± 


0.001 


435.900 




0.156 


16 


2 


912.080 


± 


0.002 


436.331 




0.213 


17 


2 


931.609 


± 


0.002 


435.855 


± 


0.180 


18 


2 


950.625 


± 


0.002 


436.627 


± 


0.146 


19 


2 


969.222 


± 


0.002 


438.012 




0.149 


27 


1 


856.964 


± 


0.002 


437.741 




0.123 


35 


1 


954.940 


± 


0.002 


440.118 


± 


0.091 



quencies, splittings, asymmetries, lifetimes, and heights 
were measured. These low-frequency p modes contribute 
to improve our resolution throughout the solar interior 
since they sample a large range of penetration depths. 
Moreover, because these modes have lower upper turning 
points in the outer part of the Sun, they are less sensi- 
tive to the turbulence and magnetic fields in the outer 
layers, which should make them extremely valuable for 
the study of the physical processes responsible for the 
oscillation excitation and damping by the turbulent con- 
vection. 

We would hke to recall that iSchoul ()1992! )'s peak- 
finding approach consists in fitting the individual-m 
spectra simultaneously by using a model in which the 



shift coefficients are introduced, while in the present tech- 
nique, the best shifts are determined first, based on the 
calculation of figure-of-merits fSec. lS.ll and Appendix [X)) . 
and then the rotation-corrected, m-averaged spectrum is 
fitted (Sec. [321). 

The development of the m-averaged spectrum tech- 
nique towards both higher frequencies and larger angular 
degrees is one of the next step to be addressed, as also 
the analysis of shorter datasets, such as the canonical 
108- and 72-day time series. 
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This work utilizes data obtained by the Global Oscil- 
lation Network Group (GONG) program, managed by 
the National Solar Observatory, which is operated by 
AURA, Inc. under a cooperative agreement with the 
National Science Foundation. The data were acquired by 
instruments operated by the Big Bear Solar Observatory, 
High Altitude Observatory, Learmonth Solar Observa- 
tory, Udaipur Solar Observatory, Instituto de Astrofi'sica 
de Canarias, and Cerro Tololo Interamerican Observa- 
tory. The GOLF and MDI instruments onboard SOHO 
are cooperative efforts to whom we are indebted. SOHO 
is a project of international collaboration between ESA 
and NASA. BiSON is funded by the Science Technology 
and Facilities Council (STFC). We thank the members 
of the BiSON team, and colleagues at the host institutes 



at each of the BiSON sites. The authors are particularly 
grateful to S. G. Korzennik for providing us with the 
2088-day MDI dataset, and to J. Schou for the MDI leak- 
age matrix. The authors thank R. A Garcia, P. Boumier, 
and W. J. Chaplin for providing estimates of GOLF and 
BiSON mode frequencies and ai rotational splittings ob- 
served in decade-long time series. The authors also thank 
S. J. Jimenez-Reyes and J. Schou for their useful com- 
ments on the manuscript, and S. G. Korzennik for help- 
ful discussions during the various stages of this work. 
D. S. acknowledges the support of the NASA SEC GIF 
grant NAG5-11703. This work has been partially funded 
by the grant PNAyA2007-62650 of the Spanish National 
Research Plan. 



APPENDIX 

FIGURES-OF-MERIT AND DETERMINATION OF THE A-COEFFICIENTS 

The best estimates of the splitting a-coefficients are obtained by maximizing the likelihood of the m-averagcd 
spectrum (see Sec. 13. ip through the calculation of a figure-of- merit (FOM). However, as shown on Fig. [1] other criteria 
to define a FOM can be used, such as the narrowest peak (i.e., the minimum mode linewidth), or the minimum entropy 
of the resulting m-averaged spectrum. In order to compare the actual mode parameters and associated uncertainties 
obtained with two different definitions of the FOM, we applied the m-averaged spectrum technique to the 3960-day 
GONG dataset by using both the maximum likelihood and the narrowest peak in the m-averaged spectrum as FOM. 
Figure lAlGl shows the corresponding central frequencies, and the odd ai, 03, and 05 splitting coefficients of the common, 
measured low-frequency p modes. The associated formal uncertainties are also represented. The two FOMs return 
consistent mode parameters within the error bars, the difference between the two being within the 3cr limit for all of 
the mode parameters. 

IMPACT OF THE FITTING MODEL USED 

As a test of the dependence of the measured frequencies on the fitting model used to describe the m-averaged spectra, 
we fitted the m-averaged spectra using three different models: an asymmetric Lorenztian profile (Eq. [2|) including the 
closest Sm zt 2 spatial leaks (hereafter called A2 and used as the reference model); a symmetric Lorenztian profile 
including the Sm ± 2 spatial leaks (hereafter S2); and an asymmetric Lorenztian profile (Eq. [2]) but omitting the 
neighbouring 6m ± 2 spatial leaks (hereafter A). Figure IB 171 shows the differences as a function of frequency in 
the 2088-day GONG low frequencies estimated using the m-averaged technique between S2 and A 2 (red dots), and 
between A and A2 (black plus signs), in both cases A 2 being the reference model. Ignoring the peak asymmetry in 
the fitting model leads to a systematic underestimation of the mode frequency as the frequency increases, the effect 
becoming particularly large above w 1400 /iHz (red dots). The differences become much larger than 3cr, for example, 
at w 1800 /iHz, the fitted frequencies between 5'^ and A 2 are abou t 20cr apart. These results obtained for modes below 
2000 /iHz confirm previous observations, e.g. iThierv et al.l (j2000[ ) who analyzed low-degree modes above 2000 /iHz in 
805 days of GOLF data. 

On the other hand, omitting the spatial leaks has no effect below « 1600 //Hz, as they become well separated 
from the main peak because the corresponding mode linewidths are much smaller than their frequency separation. 
As the frequency increases, the mode linewidths increase, and ignoring the spatial leaks in the fitting model of the 
m-averaged spectrum between about 1600 and 2000 /iHz leads to an underestimation of the target mode frequency, 
the maximum difference occuring around 1800 //Hz. The frequency separation between the target mode and the m- 
Icaks then becomes comparable to their linewidths and the lines blend together in the m-averaged spectrum. Above 
2000 //Hz, this underestimation seems to vanish. Indeed, at that frequency range, the mode linewidths are much larger 
than the frequency separation, and the first spatial leaks (at least) are totally blended into the target mode in the 
m-averaged spectrum, having a much lower impact on the frequency determination. However, the effect of ignoring 
the peak asymmetry is much larger than that from ignoring the m leaks even in the frequency range where the m 
leaks have the strongest impact. For instance, at 1800 /tHz, the effect on the frequency underestimation by ignoring 
the mode asymmetry is about seven times larger than by ignoring the m leaks. 

As an example of the other mode parameters, the right panel of Fig. IB17I shows the differences in the extracted 
mode linewidths between the different fitting models. The color code is the same as for the differences in frequency 
represented on the left panel of Fig. IB 171 Ignoring the presence of the m leaks in the fitting model leads to a 35% 
overestimation at most of the extracted linewidths in the low-frequency range showing a maximum mismatch around 
1900 /iHz. Interestingly, if the m leaks are omitted, the linewidths are underestimated below w 1600 /iHz, showing a 
maximum 10% underestimation around 1500 //Hz. On the other hand, ignoring the peak asymmetry has a very small 
influence on the fitted linewidths in the low-frequency range. However, above « 1800 //Hz, the linewidths extracted 
using an asymmetric profile are systematically larger than the ones returned using a symmetric profile. 
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Fig. A16. — Left column - l — v diagram, and ai, as, and 05 splitting coefficients for the common, low-frequency p modes measured in the 
3960-day GONG dataset obtained by using the maximum likelihood (circles) an d the narrowest peak (dots) of the m-averaged spectrum 
as figure-of-merit to determine the best estimates of the a-coefficients (see Sec. 13.11 and Fig. [TJ. Right column - Associated Icr formal 
uncertainties (nHz) of the mode central frequencies, and of the ai, as, and 05 splitting coefficients. 



DERIVATION OF ERRORS FROM AN A/-AVERAGED SPECTRUM 

The derivation of the errors of the mode central frequencies and of the a-cocfficicnts measured from the m-averaged 
spectrum technique is detailed here. 

Approximation 0/ the statistics of the m-averaged spectrum 

The m-averaged spectrum is obtained from the summation of 2Z + 1 spectra assumed to be with with 2 d.o.f 
statistics each having a different mean or signal-to-noise ratio. All of the individual 7Ti-spectra are independent from 
each other. The solar background noise is assumed to depend on m with a polynomial with only even terms (0, 2, 
etc...). The amplitude of the modes is assumed to depend on m with a different polynomial also with even terms 
(0, 2, etc.). In a first step, the Oi-coefficients are calculated to maximize the likelihood of the resulting m-averaged 
spectrum; 

Using lAppourchauxl (j2003f ). we can derive an approximation of the statistics of the summation of the 2^ + 1 spectra. 
The statistics of the m-averaged spectrum S can be approximated by a Gamma law given by: 



(CI) 
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Fig. B17. — Effect of asymmetry and spatial leaks on the fits. Left panel - Differences (in nHz) in the 2088-day GONG frequencies 
estimated using different models to describe the m-averaged spectrum: symmetric Lorentzian profile minus asymmetric Lorentzian profile 
(Eq.O both models including the first Sm ± 2 spatial leaks, i.e., S2 minus A2 (red dots); both asymmetric Lorentzian profiles (Eq. [2)l but 
ignoring the closest Sm ± 2 spatial leaks for one of them, i.e., A minus A2 (black plus signs). The differences in the extracted frequencies 
using different fitting profiles are represented as a function of frequency. Flight panel - Same as the left panel, but for the differences (in 
/iHz) in the fitted mode linewidths. 



The mean and a are given by; 



E[S]^j and = ^ (C2) 



A and vi are then derived from the mean and a as; 



A^r^and^i^:^ (C3) 



In our case, the mean E[S] and a are given by: 



m—l 



\ 



m—l 



E[S] = E /™ ^'^'i 

m— — l 

where fm is the power spectrum for azimuthal order m which can expressed as; 

uq, flj) = B„i{iy) + A,n{y, I'o, ai) (C5) 

where v is the frequency, the central frequency . a; arc the usual RitzwoUcr-Lavcly coefficients, Bm is the background 
noise, Am is the profile of the mode (the linewidth and amplitude have been omitted for simplifying the notation). 
We can write the noise as; 

Bm{y)^B{y){l+gB{m)) (C6) 

where the m-depcndence is assumed to be independent of frequency. gsiTn) is such that; 

E[Bm{i^)] = (2l + l)B{iy) (C7) 

If the correction of the has been done properly, to the first order the m-averaged spectrum is independent of the a^. 
We can write the mode amplitude as; 

Am{iy) ^ A{i^)il + hA{m)) (C8) 
where the m-depcndcncc is assumed to be independent of frequency. hA{m) is such that 

E[Am{v)] = (21 + l)A{i^) (C9) 

Then we find; 

E[S] = {21 + l){A{iy) + Biiy)) (CIO) 

and 

= {21 + 1) [A^{v) (1 + a) + B'^{iy) (1 + /?) + 2AB{v) (1 + p)] (Cll) 

with 



ra—l 

1 

a 



ra—L 

— E '^aM (C12) 
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and 



m—l 

^^^T^T^aUm) (C13) 

7n— — l 
^ m—l 

P= J2 ^A{m)gB{m) (CM) 

■m— — l 

we finally get for A and the following: 

^ A-)+B{i^) 

and vi as: 

' A^i^){l + a) + B^{i^)il + /3) + 2A{i^)B{i^){l + p) ^ ' 

After simplification we get: 

. A{i^)+Bi.) 

iA{iy)+Biiy))^ + aA^i^) + /3B^ii^) + 2pA{iy)B{iy) ^ ' 

and 

„^ {2l + l){A{u) + B{v)f 

Using the dependence observed in the GONG data, we have a « 0.17, /? « 0.035 and p « 0.075. They are sufficiently 
small such that we have: 

A«:!---^— - (C19) 

A{v)+B{v) ^ ' 

and 

w (21 + 1) (C20) 

then wc find the following statistics for the m-averaged spectrum: 

1 92; g 

P^''' r(2/ + l) (y^H+6(z.))2'+i^ ^^^^^ 
After a change of variable u = S/ {21 + 1), we have 



When we use MLE, we minimize the following 

C{i^, vo, ai) = - lnp{u) = -{21 + 1) 



ln{A{:^) + B{v)) 



(C23) 



{A{ u)+B{v))_ 

which shows that using the MLE applied to a with 2 d.o.f as prescribed by lAppourchauxl (|2003[ ) is in the case of 
the m-averaged spectrum a good approximation. It is not an approximation when averaging several power spectra 
of identical mean (or variance), i.e. when a = (3 = p = Q. Note that what we minimize is the sum over a range of 
frequency that can be approximated as: 

L{uo,ai,)^ I C{u,uo,a^)Av (C24) 



Error bars on the central frequencies 
Error bars for frequency arc derived from the inverse of the Hessian (second derivative of L) as: 

< = S- (C25) 

UVq 

iToutain fc Appourchauxl ()1994D showed that we could express the error bars as a function of the mode profile V 
{=A + B) as: 

<-^^''^'^jvh){tS'^^ ^^^^^ 

where T is t he observation time. The 2 1 + 1 factor is due to the fact that the likelihood is 2/ + 1 times larger than the 
hkelihood of iToutain fc Appourchauxl (|1994[ ) (cf Eq. 23). Eq. ((C26| shows that the error bars on the frequencies in 
the TO-averaged spectrum will be \/2l + 1 smaller than for the mean of the individual modes. In deriving Eq. (|C26p . 
we assumed that {u) =7^. This is an approximation good enough for getting the error bars on the frequency but not 
on the a,-. 



Low-frequency solar p modes in spatially-resolved helioseismic data 



17 



Error bars on the a- coefficients 
The error bars on the a-coefficients are derived from the inverse of the Hessian (second derivative of L) as: 

h 



(C27) 



As shown bv lToutain fc Appourchaiu3 (|1994[ ). these coefficients can be related to the mode profile as using Eq. (|C23P : 

1 dV^dV„., ^^23^ 



using the following property: 



dVm dVr, 



■IPi.rnim/l), 



where Pi^i''^/^) ^^'^ Ritzwollcr-Lavely polynomials. And wc finally get: 



(C29) 



(C30) 



We recog nize the error ba rs of the frequency for the m spectrum depending on the inverse of the signal-to-noise ratio 
Pm (as in lLibbrechtlll992f) : 



1 



diy. 



Finally Eq. (|C30[) can be written as: 



(C31) 



(C32) 



The errors scale like I 2 (as in lVeitzer et al.|[l993[) . If the SNR is the same for all m, then wc have cr^^^ = (2i+l)(T^^. 
Thus, by simply using the orthogonality property of the P;*^ polynomials, and as given in Sec. 13.31 (Eg. [S]). we obtain 
the following expression to calculate the error bars of the a-coefhcients in the m-averaged spectrum: 



21 



(C33) 



All terms off of the diagonal are zero. Of course, when the SNR varies with m, the off-diagonal terms arc non-zero 
and correlations appear. 
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